For practical opc applications , csplat is not suitable for intensity calculation on huge amount of sparse aerial points , hence a fast algorithm of sparse aerial image intensity calculation for lithography simulation is derived in this paper based on the theory of principal waves